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We study the dynamics of "finger" formation in Laplacian 
growth without surface tension in a channel geometry (the 
Saffman- Taylor problem). Carefully determining the role of 
boundary geometry (resulting in reflection symmetry, not pe- 
riodicity), we construct field equations of motion, these cen- 
tral to the analytic power we can here exercise. Initially we 
consider an explicit analytic class of maps to the physical 
space, a basis of solutions for infinite fluid in an infinitely 
long channel, characterized by meromorphic derivatives. We 
dynamically verify that these maps never lose analyticity in 
the course of temporal evolution, thus justifying the underly- 
ing machinery. However, the great bulk of these solutions can 
lose conformality in time, this the circumstance of finite-time 
singularities. By considerations of the nature of the analyt- 
icity of all these solutions, we do, however, show that those 
free of such singularities inevitably result in a single asymp- 
totic "finger". This is purely nonlinear behavior: the very 
early "finger" actually already has a waist, this having sig- 
nalled the end of any linear regime. The single "finger" has 
nevertheless an arbitrary width determined by initial condi- 
tions. This is in contradiction with the experimental results 
that indicate selection of a finger of width 1/2. In the last 
part of this paper we motivate that such a solution can be 
determined by the boundary conditions when the fiuid is fi- 
nite. This is a strong signal that finiteness is determinative 
of pattern selection. 



1. INTRODUCTION 

The aim of this paper is to offer some careful analy- 
sis of the initial value problem associated with Laplacian 
growth without surface tension in an infinite channel ge- 
ometry. It is well known for quite some time that this 
problem possesses asymptotic solutions in the form of a 
"finger" whose width A may be taken apparently arbi- 
trarily Q. In experiments which are performed in a fi- 
nite channel and in which the effective surface tension 
is very small a sing le width (A = 1/2) is selected |],|. 
Previous work explained how a unique steady solution is 
selected by the existence of surface tension That 
work never explored the question how this solution arises 
from arbitrary initial conditions. 

Here we consider the initial value problem, and ex- 
plain how a very wide class of initial conditions ends up 
with an asymptotic solution of an interface in the form 



of a single finger. Our initial conditions may include very 
complicated interfaces, indeed having many features. It 
is important to understand the generic dynamical mech- 
anism that is responsible for the generation of a single 
finger that becomes an asymptotic solution. We show 
that in the infinite channel the width of the finger is not 
selected (notwithstanding some recent claims to the con- 
trary§). 

We next remark that the experiments in which selec- 
tion is observed include, in addition to surface tension, 
also a second boundary condition. One can set up the ex- 
periment in various ways, like ending the material fluid 
at atmospheric pressure, or connect it to a pump that 
sustains the flow. We discuss the various boundary con- 
ditions that can be imposed on the second boundary, and 
show that some natural choices impose a unique solution 
with width 1/2. This discussion serves as an introduction 
to a substantial treatment of finiteness in a forthcoming 
paper H 

The structure of the paper is as follows: in Sect. 2 we 
recall the analytic theory of Laplacian growth in chan- 
nel geometry. The novel aspect of our treatment is the 
introduction of the equation of motion for the conformal 
map that explicitly takes into account reflection symme- 
try of the solutions in the complex plane. The significant 
virtue of our method is that it totally obviates any need 
for Hilbert transforms and so forth that are required in 
the usual, pure boundary treatment. .In Sect. 3 we show 
that analytic initial data remain analytic for all available 
time. On the contrary, conformality can be lost in finite 
time, and in Sect. 5 we discuss how this may happen in 
generic situations, and discuss how a regulation, such as 
surface tension, accomplishes the cure of this disease. In 
Sect. 4 we come to understand that any generic initial con- 
dition has a single Saffman- Taylor finger asymptotics. In 
Sect. 6 we discuss finite boundary conditions. We observe 
that the imposition of a downstream boundary condition 
may uniquely select a solution of width 1/2. In Sect. 7 we 
offer a summary and a discussion. 



ANALYTIC SOLUTIONS OF LAPLACIAN 
GROWTH IN CHANNEL GEOMETRY 
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A. The physical problem and the mathematical 
formulation 

We are interested in the so-called Saffman- Taylor prob- 
lem of determining the motion of an interface "/{t) be- 
tween two fluids of different viscosities in a Hele-Shaw 
cell. We can think of air displacing oil as a standard 
example [|j . The cell is a channel made of two long rect- 
angular plates displaced by a small distance b. We chose 
X to denote the long lateral coordinate, whereas y de- 
notes the transversal direction of the cell, < y < tt 
in suitable units, (see Fig. la). When the gap b is con- 
siderably smaller than the lateral width of the cell, and 
non-slip boundary conditions are taken at the upper and 
lower plates, then the velocity field v in the driven fluid 
satisfies Darcy's law 



(1) 



where p is the pressure field and /i the viscosity. Because 
of the assumed very small viscosity of the driving fluid, 
its pressure is almost constant (taken to be zero), while in 
the driven fluid, by virtue of incompressibility, ^ -v — 0, 
the pressure is harmonic: 



Ap = 



(2) 



The boundary conditions on the interface are determined 
by first requiring the non-penetrability of the two fluids 
in contact. This means an equality of the normal velocity 
of the interface Vi ■ n and of the normal velocity of the 
fluid at the interface. Secondly, the pressure of the fluid 
at the interface is given by the relation 



Pa 



p = an 



(3) 



where a and k are interfacial surface tension and cur- 
vature respectively. On the lateral walls {y = 0, tt) the 
normal velocity of the fluid vanishes, and as a; — *■ oo, far 
ahead of the interface, the flow is taken as uniform, par- 
allel to the X axis, and of magnitude V. Here we take 
a = 0. The boundary condition (|^) now simplifies to 
taking constant pressure on the interface: 



p — Sit the interface 7(t); a = 



(4) 



In consequence of Eq. (|^) we also need to require that the 
normal derivative of the pressure vanishes at the lateral 
walls. This means that the interface must be normal to 
the two lateral walls. As the flow is approximately two- 
dimensional and obeys Laplace's equation, in the usual 
fashion one produces an analytic function h'{z), 



h'{z) = -dxP + idyp = Vx - ii 



(5) 



where z = x + iy. (With Ap = 0, h'{z) satisfles the 
Cauchy-Riemann conditions.) The integral of h'{z) is 



h{z) = -p{z) + is{z) 



(6) 



with s{z) the harmonic function conjugate to —p{z). The 
function s(z) is typically multivalued. Consider now an 
arbitrary curve 7, and denote by $7(i) the time depen- 
dent flux that crosses 7. It is convenient to consider 
V^(i), the mean channel velocity, by dividing $ by the 
constant channel width tt. With fi denoting the right- 
handed normal 



V,{t) 



Vndi 



(7) 



where d£ is an arc-length differential. Denoting u the 
unit vector orthogonal to the plane we can write 



Vndi = V ■ {d£ X u) = V X di- u 
= Vxdy — Vydx = Imvdz . 

We observe that 

1 i I i 1 

— :vdz = V ■ d£ H — Vndi = —dp H — Vndi , 

Tri TT TT TT TT 



and 



1 /■ _ i 
Vy{t) — — : / vdz S^p , 



(8) 



(9) 



(10) 



where 5jp represents the pressure difference between the 
end points of 7. In particular, the flux crossing a curve 
of constant pressure p = Po satisfies 



1 



vdz 



(11) 



We will later need to consider sinks at finite distances. 
In preparation for this, consider two curves, 71 and 72, 
each connecting the two boundaries, and observe that 



1 

TTi 



vdz 



(12) 



where the closed curve is the boundary of the domain 
which is delineated by 71 , 72 and the two walls. If v is 
analytic in Q, V^-^ = V^^. On the other hand if there are 
sinks in 17, we can write 



ai/2 



14. 



— ' z — : 
^72 = XI °* 



(13) 



and each such pole sinks worth of the flux crossing 71 
into r2. Returning to the flux across a constant pressure 
line we note from Eqs.(pT|) and (ph that Yp„ = (5p„s/7r, or 



(14) 



Since we are assuming a flux of fluid incident from the 
left, it is evident that there must be a sink of equiva- 
lent strength either at infinity or at some finite location. 
In particular the pressure approaches —00 at the sink. 
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Provided there is just one sink, then for each value of p 
there will generally be just one corresponding physical 
curve. Then, at each instant of time, the flux crossing 
every curve of constant pressure is identical. Defining a 
complex variable 



c 



-p 



IS 



Vit) 



(15) 



where V{t) is the common value of Vp^ for all constant 
pressure lines, we now see by Eq.(p^ that h/V{t) maps 
the physical channel into an identical region of C space. 
That is, the map 

Q = h{z,t) = h{z,t)/V{t) (16) 

can conformally map the strip to itself. In particular 



h' 



V{t) 



(17) 



replaces Eq.(||). Laplace's equation for the pressure is 
written as AReC = 0, and is automatically still obeyed 
since V is purely a function of time. However, we must 
insist now that the interface is p = and not any other 
constant, not to have ReC time dependent on the inter- 
face. (The physical pressure can be any po(i)- Since only 
differences of pressure matter, we simply subtract po{t) 
to normalize to p = at the interface.) 

The relations between the physical channel, the math- 
ematical strip and h{z) are summarized graphically in 
Fig.(0)a. 



(a) 




ReC=0 



f(Q 



(b) 
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B. The analytic map / and its equation of motion 

Having defined the map h we notice that it is incon- 
venient that the boundary of its domain ^{t) is at the 
moment unknown and potentially complicated. However 
its image in C, space is elementary: ReC = 0. Accordingly 
it is natural to invert the discussion and consider a map 
/ from C to z. Assuming h is conformal {h' ^ 0) on 
the physical domain), exists, which is precisely the 
desired analytic (in the physical domain) f (f = h~^). 
From this point onwards when we say "analyticity" , we 
mean the analyticity of / (equivalent to the conformal- 
ity of h). As well, when we say "conformality" we mean 
the conformality of / (equivalent to the analyticity of h). 
Should /'s analyticity fail, then the setup is inapplica- 
ble; should /'s conformality fail, then Ap = no longer 
holds, and the dynamics has created singularities. We 
will show in Sec. 3 that / remains analytic for all times. 
The serious issue of /'s conformality is beyond our full 
understanding, but we shall illuminate the issue. f{(,t) 
describes a flow with boundary conditions p = provided 
that the interface 'y{t) develops to the interface 7(t') for 
any later time t' under transport hy v ^h' . This require- 
ment leads to the equations of motion pO|-p^ . Consider 
any point £ ^- Co serves as a Lagrangian label for the 
fluid point at /(CoiO). However, /(CoiO not this same 
fluid point at time t. Define 



z = Z{(o,t) 



(18) 



as the moving point labeled by Coj with Z{(q,0) — 
/(Co,0). There exist at each time t a map Co ^ C(Coj*) 
such that Z(Co,t) = ,f(C{(Q,t),t). By definition 



ZtiCo,t) = f'Ct + ft^v,+n 



V{t)h'{f{C,t)^ 



V(t) 



net) 

Accordingly ( satisfies the differential equation 

l/fCt + /*/' ^ Vit) . 



(19) 



(20) 



We observe now that the boundary j{t) flows into itself. 
Thus for each Co on the boundary (ReC = a.t t — 0), 
C(Co, t) is also on the boundary, and so ReC(Co, ^) = for 
ReCo — 0. But then ReCt — and taking the real part of 



j) we obtain the usual result 1 10 

Rc(/t/') = v{t) 



on ReC = 



(21) 



In previous work V(t) is taken identically equal to 1. 
Eq. (|l]) generalizes the standard result to the case of vari- 
able flux which is necessary with finite boundary condi- 
tions. 



FIG. 1. (a) The physical channel in z space with the in- 
terface satisfying no-flux on the lateral walls is mapped onto 
the mathematical channel in C = h{z) ~ —p{z) + is{z). (b) 
Using the reflection symmetry the problem is "doubled", the 
interface in the doubled channel satisfies periodic boundary 
conditions, and the expected physical solutions is two fingers. 



C. Exponentiation and reflection 

Geometrically, the strip SI, ReC > 0, < ImC < tt, 
is less felicitous then a circular domain, with "infinity" 
(ReC ^ oo) the point at infinity. It is natural to write 
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w 



giu), /(C)=ln.g(eO 



(22) 



so that is the entire upper half-plane minus the unit 
disc (|u| > 1, Imu > 0), see the shaded region in Fig.(||). 
Under /'s exponential conjugate, g, maps to Imw > 0: 




FIG. 2. The exponentiated domain and its mapping g{u) 
to the exponentiated physical strip 

g's boundary value on |Reu| > 1, Imu — > 0^, is real, 
and just the exponential of pressure along the walls. 
With no shocks, p{x,0) is continuous. So g takes real 
rays continuously to real rays. As such it immediately 
defines its analytic continuation g{u) = g{u) for Imu < 0, 
and this reflection symmetric continuation is analytic in 
the entire u- plane for |u| > 1. With In and exp both also 
reflection-symmetric, so too must be /: 



/(C) = /(C) , 



(23) 



as well as /' and ft. We shall have numerous occasions 
to utilize this symmetry. Note that the continued / has a 
domain equivalent to a doubling of the physical problem 
such that the interface and fluid now enjoy analytically 
periodic boundary conditions from — tt to tt. Of course, 
the physical problem remains in the half strip, but this 
extra symmetry dictates profound restrictions on /'s an- 
alyticity structure |l4[. 



D. Reflection symmetric equations of motion 

Let us now utilize reflection symmetry to transform 
Eq.(pTl) into a field equation valid throughout the body 
of the fluid. Using (|2^) we can rewrite Eq.([2l|) in the 
form 

2U(<) = /t(C)/'(C) + /t(C)/'(C) , (24) 
and noticing that Rc^ = implies that C, = — C we have 

/t(C)/'(-C) + M-Of'iC) = 2U(t) . (25) 

It is important to note that this equation holds not just 
for ReC — 0, but for all ( for which f is analytic. The rea- 
son is as follows: solve for /(C, t) from Eq.(pi|). Compare 



to the solution of Eq.(p5|), which we denote temporar- 
ily as f{C,t)- The two functions f{C,t) and f{C,t) are 
the same on the boundary Re^ = 0, and since they are 
analytic they are the same functions. We can thus use 
(H) instead of (|l]) for any C- Eq.(|||) relates /(C) to 
/(— C) which will prove to be of significant value for the 
forthcoming analysis. In particular, when / has the el- 
ementary analytic behavior /(C) C as ReC — > cxd, it 
follows from ( p5| ) that / is analytic as ReC — oo. (This 
is the usual case for the consideration of a sink at oo.) 



E. A class of solutions 



There 



exist, sprinkled throughout the literature 
3|jl5[, considerations of a class of solutions to (|25|) 



with a sink at oo, with moving singularities, most natu- 
rally written in u-space as 



5(w) 



e^(*)u 



n 



1 



afe(i) 



(26) 



It is easy to see that for such g the corresponding ft and 
/' simply have poles in e^. From the form of (25) we 
then see that the class of functions g exhausts the so- 
lutions of ( p5| ) within the class of rational functions of 
e^. We will verify (as has been commented upon in some 
of the literature [0|) that the equations of motion are 
satisfied with ak constant in time. Contrary to some of 
the literature we consider the ak arbitrary complex num- 
bers rather than just real integers [p^ . Further, paying 
attention to the analyticity structure imposed by refiec- 
tion symmetry, it follows that /3(t) is real, \ak{t)\ < 1, 
and for each complex {ak{t), ak} pair, there must also 
be a corresponding pair with both conjugated. Those 
g^s with just one a (which is real and constant in time) 
comprise the Saffman- Taylor solutions. Numerical explo- 
rations in the literature appear to indicate [p^ that richer 
5's replicate a large class of interface motions. (From a 
numerical point of view a rich enough collection of a's 
surely serves as a basis.) Some of our comments will fol- 
low directly from the equations of motion, whereas others 
pertain only to arbitrary solutions of the form ( p6| ) . By 
the conjugacy (22) we now have for / 



N 



_ pCkit)-C 



(27) 



fe=i 



where afe(t) = exp(Cfe(t)), Re(k{t) < 0. We shall s usp end 
until Sect. 3 the equations of motion implied by (g5|) for 
the parameters of /, and further comments on the origin 
of class (|7|) in subsection G to follow. 



F. Elementary flow solutions 

Our equations of motion allow for the ready produc- 
tion of a variety of solutions. Consider first the class of 
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solutions with a fixed shape interface uniformly translat- 
ing in time, so that we can take V{t) = 1. The interface 
is 



Zint{s,t) = f{is,t)=Pit)+F(is) 



(28) 



Using (M 



(29) 



But then, = A, a constant, and f3 — t/X. From (|29| ) 
we deduce 

F(C) - F(-C) = 2AC , (30) 

or /(C, t) - /(-C, t) = 2AC. Setting C = is, we find 

/(^s,^)-/(-^s,<) = 2^As. (31) 

By reflection symmetry f{~is) — f{is), i.e. Im/(js,i) = 
As, or 



yis,t) = As . 



(32) 



This solution represents an interface that occupies a con- 
stant channel fraction. A, sensible for A € (0,1]. This 
also implies that the interface is a graph of {x{y), y), and 
so conformal for those A's since |/'| > \y'\ = A. Writing 
F = AC + S(C), Eq.(|0|) is 



(33) 



E{-0=E{C) . 
The general solution is 

f = \ + K + E(C) 



(34) 



with E an even function of Q compatible with full channel 
width: 

X + — {E{-p + iT:)-E{-p))^l, (35) 
111 

for those curves of constant pressure within the physical 
fluid [p < 0) which span the free channel. Since E is even, 
if E has a singularity at (^k with ReCfc < 0, then it must 
also possess an identical singularity at Re^fc > 0, and so / 
is no longer analytic in the entire physical fluid, implying 
the existence of stagnation points or inflnite velocities 
within the fluid itself. Within the class of solutions g, 



/ = | + AC + ^afcln(l-e«'=-^)(l 
ReCfc > (by definition) . 



(36) 



It is expedient, to facilitate the generation of elemen- 
tary solutions, to consider those solutions which are "pe- 
riodic" over the physical channel: 



/(C + z^) = /(C) + «^ 



(37) 



Defining now double-width variables 

C = 2C — in] z = 2z — iiT 



(38) 



so that the physical channel <Imz < is mapped to 
— TT <Imz < TT, we consider 



^ = 2/ ^ -.../(C) 



(39) 



/ now has 27ri periodicity, and is again reflection sym- 
metric since / is: 

/(C)=2/fi±^V- = 2/f^V- (40) 



-2/ ^-.. +.. = 2/ ^ /(C) 



Conversely, any 27ri periodic reflection symmetric / de- 
termines a iri periodic reflection-symmetric /(C)- (Such 
an / is too symmetric in that the upper and lower chan- 
nel walls have perfectly synchronized flows along them, a 
condition not enforced by the physical data.) The virtue 
of this device is that an / with just one symmetric pair 
of singularities produces an / with two symmetric pairs 
of singularities, and so embraces more complicated flows 
that does just one symmetric pair. 

So far as time-dependence is concerned it is natural to 
take 



/(C,t) = i/(2C 



lTT,2t)+i- 



{i=2t) , 



so that 



f'iC,t) = fiC,i) , MC,t)^Mc,i) 



(41) 



(42) 



Also, by 27ri periodicity of /, f'{~C,t) = /'{-(, i) and 
similarly for ft- With V{i) = V{t), the equations of mo- 
tion are covariant. So, all we need do is forget the tildes, 
solve the simpler problem, and then conjugate it back to 
the physical space. In the sequel we neglect the last step 
as an ellipsis the reader readily can fill in. There are rea- 
sons to be wary of the extra symmetry in the solutions, 
and we shall point out this fact when it arises. 
The simplest possibility of (pq) is one Cfc = ""i: 



/ = - + AC + a ln(l + e-'^)(l + e'^) 
A 

= ^ + iX + a)C + 2aln(l + e^^) . 
A 



(43) 



For full channel width for ReC > 0, we require X + a = 1, 



/ = - + C + 2(1 - A) ln(l + e-'^) (/'(-oo) = 2A - 1) 



(44) 
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These are precisely the Saff man- Taylor solutions for a 
finger of width A, and the only such solutions of our class 
analytic for all ReC > 0. 

The next simplest solution is 



/ = ^+C + 2(l-A)ln(l + e-^) 



aln(l 



)(1 



=-c-+c 



) (C*>0) 



(45) 



Notice that for — C* < Re^ < C*: the channel is fully oc- 
cupied with flux, but for Re^ > C* the entire asymptotic 
flux occupies only the fraction 1 — a of the whole channel. 
The most interesting such case is a = 1, when all flux is 
sunk in the point /(+oo) = Zsink^ 



Zsini.it) = J , /'( + 00) = . 



(46) 



The constant pressure contours for > p > — C* span 
the channel, but those for p < — are closed curves 
surrounding the sink at Zgink- The separatrix between 
them at C = C* + is has the form 



t 
A 

ln(l 



+ C* + is + 2(1 - A) ln(l + e-^*-*") 

-2C,-ts-^ _ ^^^^^ ^ts^ _ (-47) 



For C* > 1 



-+ + is -h-i{l + e'') , 
A 



(48) 



and is an upstream-pointing Saffman- Taylor finger of 
width 2(1 — A) = 1, or A = 1/2, by comparison to (44), 
see Fig.@. 




This is amusing, and suggests that with an enforced 
symmetry between C and — ^, then ( |4^ ) would be 
f'{—oo) — 0, so that if (^J) were asymptotically valid, 
then a 1/2 width finger for the interface itself could be 
implied. We will return to these matters later. 

Let us now consider some elementary solutions with 
variable flux V(t). The simplest solution arises when 
the pressure profiles are a function of x but not of y, 
p — p{x,t). This implies, by analyticity, that 



/ = m + c 



(49) 



where the coefficient 1 of C specifies a channel fully occu- 
pied with flux. We now pose finite boundary conditions, 
namely that p = Pa (atmospheric pressure) on x_= L. 
We denote pg = \pa\, and write 2L = /(C) + /(C)i on 
2pg/V = C + C Accordingly 



2L = /(C) + /(|^ 



C) 



(50) 



or L = P + pg/V. Using the fact that in this case $ = V, 
we need to solve the differential equation 



(3P-Lf3 + pg . 



(51) 



The solution is — 2Lj3 + 2pgt = const. Choosing initial 
conditions such that the interface is at zero when i = 
determines Xint = (3 = L — ^ LP- — 2pgt. For the velocity 
we find 



V^I3 = 



Pg 



^lJ^2pit 



(52) 



This makes it clear that any genuine physical determi- 
nation of the flow (including the flux) requires boundary 
conditions at finite pressures, rather than at p — > — cxd. 



G. The class of solutions (27) as a resummation 

Before turning to the discussion of analyticity and con- 
formality of the solutions, it is important to reflect on the 
meaning of these solutions. It is straightforward to ver- 
ify that if we start with a Saffman- Taylor solution with 
width 1/2, that is. 



h = 2t + ln(e 



c-c 



1) 



= 1 



(53) 



FIG. 3. Schematic stream and pressures lines in a channel 
with an advancing finger (shaded region) and a sink (denoted 
by X ) . The stream and pressure lines are indicated by dashed 
and continuous curves respectively. The pressure line with 
p = p* separates pressure lines that connect the walls (like 
Pa) and those which form closed curves around the sink (like 
Pb)- The stream lines s — vr, s = —n graze the walls y = ±7r 
respectively from a; = — ootoa:: = oo where they curve around 
and return along y — 0'^ and y — 0~ respectively from x = 00 
to the sink, forming a branch cut. 



and perturb it to say fo+ef, then the equations of motion 
are satisfied with errors O(e^), for each 



n= 1,2,. 



(54) 



all analytic as ReC — » 00. That is, the basis for perturba- 
tions to /o grow exponentially in time. We can construct 
from these, "center manifold" perturbations that grow 
slower than exponential in time, by forming the infinite 



sum 
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/a = -^(-a)'7n = ln(l + ae2*(l 
1 

= ln(l + ae2*)+ln(l + e^=(*)-^) 
Cs = -C-ln(l + e-27a) . 



)) 



(55) 



(For a = 1, Cs = 0- Thus, class ( p7| ) viewed as /o plus 
fluctuations, is simply an infinite resummation of the un- 
stable manifold spanned by (|4|). But this makes it clear 
that generic functions are to obtained, reciprocally, by 
infinite summations over class (|27|), they behaving very 
differently from a finite sum. This leaves as open the 
issue if conclusions based on class ( p^ ) have any unre- 
stricted applicability. This question is beyond the scope 
of this paper, and will be addressed in a future paper. 



3. ANALYTICITY AND CONFORMALITY AS A 
FUNCTION OF TIME 



(Had we taken the ak time dependent the first sub- 
dominant asymptotic term would have revealed that in- 
deed the ak are constant.) Recognizing that this equa- 
tion reads df{—Ck,t)/dt = 0, we introduce the points Zk 
in physical space. 



Zk 



f{-Ck,t) , ife = 



(64) 



(Notice that by refiection symmetry Zk, as the Cfc, come 
in complex conjugate pairs.) With ReCfe < we see that 
— Cfe is in the interior of Q, and therefore Zk is within the 
moving fluid. Thus Zk are asymptotic stagnation points 
of the flow as the singularities approach the interface. 
This result was observed for the flrst time in More 
correctly, the fluid stagnates ak In 2 upstream from Zk as 
ReC/c ^ 0, SLt ( — ilm(k- Finally we rewrite Eq.(|6^) in 
the form 



Zk = Pit) -Ck+Yl - ^'^'^'^'') 



(65) 



A. Solutions of the equations of motion for the class 
g and asymptotic stagnation points 



Considering the solution (27) we compute 

N 



ttfcCfc 



fc=i 

N 



;C-a- - 1 



ak 



,C-Ck 



k=l 



1 



(56) 
(57) 



It is useful to consider these equations in various limits. 
Consider flrst C — > ±oo on the real axis: 



C -oo : = /3 + E 



(58) 
(59) 

(60) 



Substituting in the equation of motion ( |25[) we flnd 

(2-^afe)/3 + ^afea = 2(t-to) , (61) 

k k 

with to a constant of integration. Next consider the 
asymptotic behavior C Cfc■ 



C ~ Cfc : ft 



OtkCk 



ak 



(62) 



Substituting in the equation of motion (|25|), after divid- 
ingby/'(C)/'(-C), wefind 



a/'(-a,t) + /t(-a,t) = o 



(63) 



B. Necessary Conditions for Asymptotic Analyticity 

The discussion of analyticity is largely independent 
of the explicit form of solutions (27). Infinite channel 
asymptotics implies that 



/ = /3 + C + F(C,t) 



(66) 



with the first two terms describing a uniformly translat- 
ing fluid, and the last term its decoration, (exponentially) 
vanishing at +oo, so that 



F, F', Ft ^0 as ReC 



(67) 



Much of the discussion relies just on this. Moreover, to 
exhibit the moving singularities, we write 

F{c,t) = Y,Fk{c-Ck{t)) , Ft^-j2CkFuc-Ck{t)) . 

(68) 



with each Fk obeying (|6^), and by the analyticity of /' 
at— oo. 



Ki^oo) = -ak 



(69) 



Granted this level of detail of /'s form, we have by the 
equations of mot ion ([2^), with f ^ f3 + ( as Re^ +oo. 



/t(-oo) + /3/'(-oo) = 2 , 
(l + /'(-^))/3 + Ft(-^) = 2 , 
(l + /'(-^))/3 + E«fea = 2, 



(70) 



(71) 



(72) 
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(l + /'(-^))/3 + ^afea = 2(<-io) . (73) 

Let us write 

l + /'(-oo) = 2-^afc = 2A. (74) 

Should all work well, with no impediment against all the 
Ck a pproaching the interface Re^ = as i — > oo, then 
(|73|) implies 



,4 



(75) 



It is now straightforward to show that the physical case 
requires 



and hence by (p 



< A < 1 , 



l/'(-«))l<l 



(76) 



(77) 



To see this, let us determine C(Coi for ReC — * +oo. The 
differential equation for (, Eq.(EQ) reads here 



1 - Ct + /? ■ 



Integrating 



(78) 



(79) 



A given, far downstream, (^o lies on a line of constant 
pressure, imaged by / at t = into a curve in z-space of 
constant pressure. Should C(Co,t) at later times be fur- 
ther downstream, then by (^) and ( |67|) this same fluid 
particle lies on successively /Matter pressure curves, so that 
a flattening profile propagates upstream towards the in- 
terface. Of course, precisely the opposite must occur, 
and so by (|7|) 



1-Y<0, orO<A<l 
A 



(80) 



WithCof (f|), we now know the trajectories of far down- 
stream particles: 



/(C,t)^/3 + C~t + Co 



(81) 



just reflecting uniform unit flux of the fluid. On the other 
hand, with the interface at Re^ ~ and ReCfe — ^ 0^, 
F is bounded outside sufficiently small disks about the 
Cfc's, and so F{is,t) is bounded outside sufficiently small 
intervals of s. Observing Eq.(p6|) and realizing that for 
long times — * 0, (and see a formal proof in the next 
section), we reach the conclusion that 



"4 



(82) 



and so the interface is moving at velocity 1/A > 1. Con- 
servation of flux then imposes that the net width of the 
moving interface is just A times the full channel width, 
which to be physical (not to fail /I's conformality) must 
lie between and 1. It is thus clear that of all possible 
/'s, only those obeying ( |77| ) meet physical boundary con- 
ditions. We next show that for just these /'s analyticity 
is never lost in finite time. 



C. Absence of Violations of Analyticity 

We prove now that the solutions do not lose analyticity 
so long as ( |77| ) is obeyed. / can fail to be analytic in the 
physical region in two possible ways: either \(3\ diverges, 
leaving / defined nowhere, or the moving singularities 
cross ReC = into the physical regime. By presumption, 
/3 is finite and ReCfc < at t = 0, so that until such a 
disease occurs, ReC/t are bounded from above. We shall 
now demonstrate /'s analyticity for all future times by 
exhausting the possibilities. First, consider that some of 
the Cfc's tend to — oo at some finite or infinite time. Call 
the set of all such indices k S, and its complement, S, 
then the set of singularities remaining bounded. We show 
S is empty without proviso. According to the equations 
of motion for the singularities, for Cfc — oo, by ( |66[ ) and 



) 



Zk 



/(-a) ~ /? - a 



and produces a disease. But then 



(83) 



(84) 



and by (|73|) and (174|) 

2t (3{2 ~^ak) + f3^ak + bdd = /3(2 - ^ Uk) + bdd 
k s s 

(85) 

Should 2 — a = this can (potentially) occur at ffiiite 
time, else as t ^ oo, 2 — J^s ct < {(3 ^ — oo). That is, 
we require 



(86) 



^«>2. 

s 

Now, consider (^k with any k S. From ( |6^ and (|6£ 

(87) 

or, with (3 — > — oo 

J2 FfX-Ck - 0) +0O, for each keS . (88) 
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Eq.(p8|) can be satisfied only if Cfe + ^ 0, ±27ri so that 
becomes singular, i.e. when all the ReCfc 0~ , k € S. 
Define 



Sk = {£Kk + Q^0,±2TTt} 



(89) 



Clearly, fc G S'fc, so that no Sk is empty, and the S'^'s are 
a partition of 5. Then, (RSh is 



(90) 



Here we use a detail of of (^): the logarithm tends 
to — cxD, and so (|9(j) requires 



(91) 



Re ^ < , 



where Cfe = Cfc- Summing 
have 



Q) over distinct 5*^, we then 



y^Q; = Rey^a<0, 



(92) 



but this contradicts (pG]). Finally should S be empty, 
(all Cfc — > — oo), ( ^ is t ~ /3, impossible since t > 0. We 
thus conclude that no can go to — oo for any t > 0. 
We are now left with the circumstance that all (k are 
bounded. Should /3 be finite, then it is impossible for any 
Cfc to cross Re( — 0. This is contingent upon Fk imaging 
an arbitrarily small disc about (^k with arbitrarily large 
modulus, such as is the case for Fk of ( ^ ) with Rea^ ^ 0. 
In this circumstance, and with /3 finite, |/(Cfc)| ~* oo, 
/(~Cfc)l = kfcl =fimte, and so (^k cannot approach — C/c, 
i.e. Re(k cannot approach 0. Accordingly the discussion 
has contracted to ReCfe 0, \(3\ oo, for which ( |85| ) 
reads 



2t 



/3(2-^afc)=/3(l + /'(-(X3)) 



(93) 



Observation 1: Unless /'(— oo) = — 1, then / remains 
analytic throughout the physical region for all positive 
time, t can only continue to +oo if |/'(— oo)| < 1. 

The curious part of this result is that if |/'(— oo)| > 1, 
then / remains analytic, but the Cfc never get to cross 
RcC=0. This happens because then / must lose confor- 
mality because ft has diverged (i.e. /3 and Cfe's become 
infinite) although / is finite and analytic at such an in- 
stant. Should ft diverge, then the equations of motion 
imply that /' = at some point on the interface, i.e. that 
u ^ oo at some point of the physical fluid. What our ob- 
servation says is there must surely be finite-time singular- 
ities should |/'(— oo)| > 1. These singularities represent 
precisely an impending violation of boundary conditions. 
To see this, recall that our wall boundary conditions are 
Xs{p,s) = on s = 0,7r. By reflection symmetry, / has 
been analytically continued to — tt < s < tt, rather than 
just the physical < s < n. However such an / does not 
necessarily map s > to just y > 0. It will do so pro- 
vided that / is conformal. Should / not be a contraction 
at —oo, it is inevitable that, as the singularities approach 
the interface, / will begin to map across half channels, 
and so, such an / fails physical boundary conditions, and 
so is not an admissible solution. (To see this, by Eq.(p2[), 
g' / g — f /u, or g ^ (-°^) as u — > 0- i.e. everywhere 
within the radius of that ak{t) with smallest modulus. 
That is, with argu = <j), argw ~ /'(— oo)0, and exceeds 
half channel width unless |/'(— oo)| < 1.) Consequently 
our observation reads that all admissible solutions un- 
der physical channel boundary conditions are free from 
finite time failures of analyticity, so that the conformal 
machinery we are employing is perfectly applicable. The 
difficulty is, while |/'(— oo)| > 1 always leads to finite 
time failures in /'s conformality, the same disease can 
arise for admissible solutions as well. This difficulty is 
profound, and we will expand upon it in the next sec- 
tion. 



Should /'(— oo) = —1, then /3 can diverge at finite t, 
and generally a finite-time loss of analyticity can occur. 
Otherwise we have just t —>■ +oo , and so, 



/'(-w)>-l, /3^+oo 



(94) 



(95) 



(Since \(3\ — > oo, every Cfc crosses Re^ — "simultane- 
ously" Sit t ~ +00.) Should case (^) hold the argument 
leading to (|9^) is appropriate, save that S includes all 
the Cfe's, and so ^ Offc < 0, i.e., f'{—(x) > 1, a contra- 
diction to (|9^). Thus, only case (|93) remains, in which 
case, the same argument leadin g to (|9d) is now > 0, 

i.e. /'(-oo) < 1, which with @ is |/'(-oo)| < 1. To 
summarize what we have now demonstrated, 



4. THE ASYMPTOTICS AND THE 
EMERGENCE OF ONE FINGER 

We are ready now to establish a significant result: For 
t — > oo the physical channel supports one finger. In the 
doubled channel, that is considered in all the literature, 
this of course means two fingers with a stagnation point 
on the symmetry line. In particular, one finger in this 
geometry is a physically incorrect result. Mathematically 
this is equivalent to the following 

Observation 2: For solutions that do not lose confor- 
mality for t ^ oo, and for which all ak are generic with 
Imak ^ 0, all (k 0, ztin. 



Demonstration: Consider one of the terms in the sum 
say ak In (l — e^''^^) . Consider a circle of radius |e| 
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around the singularity, ^ = ^j. + e, where e = lele'*, |e| is 
small, and — tt < cj) < it. We investigate the image of this 
small circle under /. For small |e| 



/ ~ P{t) + Q;fe(ln |e| + i(t)) mod ± iir 



(96) 



For a given |e| this image is a line of length 27r|Q;fe| per- 
pendicular to the direction of ak (see Fig. 3). Consider 
next the punctured disk around the singularity with ra- 
dius |e|. This disk is imaged onto a series of slanted strips 
ordered along the direction of ak with the above width 
as shown in Fig.0. 



FIG. 4. The mapping of a neighborhood of one of the sin- 
gularities 

As time progresses and /? increases, this series moves 
to the right, entering the physical domain. (Remember, 
this is a neighborhood of a singularity, not a singularity!). 
Eventually, when |2ReCfc| < |e| the image of a very small 
|e|-neighborhood of Cfc must include the stagnation point 
Zfe, the fixed image of — Cfej smd must for all future time 
always include it. But as t increases, the fixed width 
strip containing Zk is moving arbitrarily far to the right. 
This would seem impossible. The only way to resolve 
this conundrum is that 



Ck (t) ^ or ± ITT a.s t 
When this happens e'^'' 



» . (97) 
and Ck (mod 27ri if neces- 



sary) is also in the |e|-neighborhood, and its term in ( |27] ) 
must be also included in (y&j. But 



aln(l 



aln(l 



) 2Realn(l - e 



(98) 



and so the strips begin to rotate to the horizontal, and 
sufficiently so that Zk always remains within its strip. 
( P?! ) is the only way to have Re^fe 0. 

It should be stressed that even though an asymptotic 
finger solution is emerging, its width is in no way selected. 
Moreover, over the duration of an actual experiment, not 
all the singularities need be yet in asymptotic proximity 
to ReC — 0. That is, a subset of (k's of the same order of 
magnitude will carry out migration toward or itt, while 
other sets with extremely large -Re^^ are still very far 
from asymptotic behavior. 



To see this, consider as usual / of (27) built from those 
Ck that are sufficiently close to one another, and consider 
an exponentially small perturbation, built from C'k with 
ReCfc < -a; < ReCfc, x > 1: 



E«feMl-e^^-'^). 



Clearly, near the interface the last term is 0{e 
ing the equations for C^, 



z'k^fi- Ck + 0{e-^ 
Writing the /3 equation 



) 



Rez[. > P + x 



(99) 
Writ- 

(100) 



{1 -J2(^k)(3 + p +J2^kCk 

+ (E"fcC^-/3E"fe)^2t-fc . (101) 



Multiplying each equation of ( |lOO[ ) by the a'^ and sum 
ming, we see, with k — '^a'^ 



(1 - E + /? + E "fcCfe = 2* + 0(e-") , (102) 
while the equations for each Qk are 

P-Ck + Y,'^t^<^-^'^''^^') + 0{e-'-') . 



Zk 



(103) 



We see then, by ( 102^ and ( |103| ) that the near set of Cfc's, 
for X large enough, behaves exactly as a full system, and 
forms a finger of w idth A = 1 — and (3 = t/X. 

But then, by ( p^ ), for t < Ax/2, ReC^ < -x/2, and the 
C'k still play an exponentially small role up to very long 
times during which a A-finger is propagating. Ultimately 
for t ^ Ax, the C'k become asymptotic, and the finger 
metamorphoses into a A' = 1 — (Y^ak + X^'^fe)/^- So, 
from a physical viewpoint, with "^fc + '^fc ~ 1' 
finger will become A = 1/2 long after the experiment is 
over. But f'{—oo) = 1 — — = 0no matter how 
large x may be. We see from this that /'(— oo) = is 
insufficient to determine A = 1/2. The width is deter- 
mined just by the CfcS near enough to ReC = 0. Indeed, 
in ( p9[ ) the X '^'k '^^n be chosen arbitrarily with no con- 
sequence during the experiment. Note that (^) is the 
transparent implementation of the manipulations of [Q. 
It is clear that |0] then has no significance for selection 
of a 1/2 finger. (By choosing the X o^' arbirtrarily, the 
argument of would then show the "selection" of any A 
whatever.) These comments pose a limitation on Obser- 
vation 2 as well: The larger the number of Cfc's we choose, 
the more semi-asymptotic regimes the initial conditions 
can be chosen to determine, so that in the limit as this 
number diverges, no final asymptotics need exist. 



VIOLATIONS OF CONFORMALITY: 
TIME SINGULARITIES 



FINITE 



The discussion that follows critically assumes that / is 
in the manifold (E^) . Should the velocity of a fluid point 
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diverge, h' ^ oo and correspondingly /' — > 0, and so the 
map /, while remaining analytic, has locally lost confor- 
mality. As we observed at the end of Sect. 3.c, should 
\f'{—oo)\ > 1, for infinite channel flow, then there must 
occur a violation of conformality. However, this circum- 
stance is trivial, in that we realized the failure occurred 
on the walls of the channel, and hence is not a solution 
under our boundary conditions. Regrettably, this is far 
from the only way in which such violations occur. How- 
ever, it is not too hard to determine when /'(C) = for 
some C within the physical fluid. Since for infinite chan- 
nel flow /'(oo) — 1, and so long as /'(C) 7^ for ReC > 0, 
then I /'(C) I > throughout the physical region. But 
then I /'I, since /' is analytic, has its minimum value over 
a region on the boundary of that region, which now evi- 
dently means the interface at ReC=0. It now follows that 
if a failure of conformality occurs, it must first appear on 
the interface, and so the condition for a finite-time sin- 
gularity is 



/' = = 1 + 



fiis) = 



(104) 



for some < s < n (but not a boundary violation at 
or 7r). 

To understand what happens, consider the behavior of 
the Saffman- Taylor solutions, that is = 1 in (p7|). By 
(1) 



/3 



2t 



apo 



2- a 



(105) 



where the one singularity, po, is real as is its correspond- 
ing a. Then, by (|65|), with zq = 0, 



0= t^+aln(l-e2''«) 



or 



t = Po - T^a(2 - a) ln(l 



(106) 



(107) 



(|107|) is soluble for pq as long as t'{pQ) = dt/dpo ^ 0. For 
Po ^ — oo, evidently t ^ po and t can increase from the 
far past. For t' = 0, generally, t will have then a maxi- 
mum, and so Po oo, inducing a finite time singularity. 
But 



t' = 0=1 



a{2 - a) 
e-2po _ 1 



-2po 



{l-af = [fi-oo)Y 



(108) 

= |/'(-oo)| . 

(109) 



We now see that with po < 0, t' = is impossible for 
/'(— oo)l < 1, but certain otherwise. That is, time 
"locks" for the inadmissible cases, and only these cases. 
With p — > oo (but f finite and analytic), ft ^ oo, and we 
seek a C on the interface where /'(C) = 0. But 



;C-PO _ 1 



aC - 



(l-a)e'' 



r(-oo) 

l/'(-oo)| 



±1 



(110) 



(111) 



This is the general nature of the failure of conformality. 
At some time, to, time locks and various of f3 and Cfe di- 
verge, and hence /( diverges, although / can be perfectly 
finite. Simply, the system (|6l|) , (|6^) becomes locally non- 
invertible for the Cfc's- 

We noticed in Sect. 2. F that Re/'(C) = A for our trans- 
lating solutions so that |/'(C)| = A > 0, and the in- 
terface is conformal, and any violation of conformality 
is the interior representation of sinks or sources. But 
Re/'(is) = y'{s), ot y = As, and so x{s) = X{y), and 
the interface is a graph of x on y. Generally, a graph 
with finite X'(y) (i.e. a differentiable graph) won't fail 
conformality. 

For the class of solutions in Sect 2.E with all > 0, 
it is easy to see that y'{s) > A, and so y is monotone in 
s, and hence the solution is a graph with |/'(C)| > A and 
so always conformal for < A < 1. The only possibility 
for a failure of conformality is with complex a's and the 
interface, a graph in the far past, about to become not a 
graph. (Indeed, the generic rotation mechanism for the 
emergence of one asymptotic finger is Imak for all k.) 
Real time singularities thus can arise when a "balloon" 
(not a graph) is about to form. 

With ak = ak + ibk, Re/' (is) is modified with bk ^ 
by terms exponentially small when the singularities are 
far from ReC=0, that is ^ e^'^'''' . Thus conformality can 
fail only when singularities enter the rotation mechanism 
of Sect. 4, which turns fingers into balloons. In particu- 
lar, this is definitely beyond the perturbative regime, and 
when the nonlinearities have become very strong. The 
usual linear stability analysis, with temporal exponents 
proportional to wavelength, simply means that fluctu- 
ations very rapidly bring the solution into the strongly 
nonlinear regime. Indeed, rather than infinitely wrinkled, 
distorted interfaces, if the rotation mechanism can work, 
a smooth single balloon is the consequence of the nonlin- 
earities, provided class ( p7| ) obtains. We will have more 
to say about this in [H , as it transpires that the unstable 
bhavior of the infinite channel is physically significantly 
wrong. Any inspection of early interface structure - say 
in Saffman- Taylor's original paper - reveals that it is bal- 
loons that are pervasive, and not graphs or fingers. Let 
us consider the simplest balloon. 

With the channel-doubling conjugacy of Sect. 2.F im- 
plicit, consider the solutions with one C^: 

f = P + C + a ln(l - e'^""'^) + aln(l - e^""'^) 
a = a + ib; a,b>0; Co = + C > ■ (112) 
The equation of motion zq = /(—Co) is 



zq = P - + a ln(l — 



a ln(l — e 



2Co^| 



zo = P + i + iV + a ln(l - e^^C) + aln(l - e-2?-2«'?) . 
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Its imaginary part is 



6. FINITENESS 



yo = + aArg(e2« - g-^"') - - In , (113) 



or 



yo — r] + a tan 



-1 



sin^ r] 



sin 2?? _ b 

e2?- cos 277 2 ^ sinh^^' 



(114) 



with the usual principal value of tan~^ correct, yo is here 
half the distance between the two "stagnati on" points 
zq and zq. The level curves of the RHS of (114) with 
< 2/0 < I" are precisely, for each j/q, the trajectory 
curve of a Co, however it be parametrized by t. There are 
three types of trajectories that connect to ^ oo (with 
y ^ 2/0 as ^ -> oo): 

(i) tan~^ f < 2/0 < TJ": the trajectory nionotonically (in 
increases from 2/0 to tt as ^ — > 0. 

(ii) atan~i f ~ I ln(l + fr) < 2/0 < tan~^ f : the trajec- 
tory moves from 2/0 to tt as ^ — > 0, initially to lower i] 
values, and with a unique minimum. 

(iii) < 2/0 < atan~^ | - 61n(l + p-): the trajectory 
monotonically flows from 2/0 to as ^ —> 00. 

Trajectories of types (i) and (ii) rotate to tt, and the 
"walls" at ±7r are closed to flow, with a balloon symmet- 
ric about 2/ = moving down the channel. Type (iii) has 
Co and Co both rotate to 77 = 0, blocking flow along y — 0, 
with fluid advancing along the ±7r walls. By Sect. 2. F, (i) 
and (ii) have blocked flow at 2/ = 0, tt with the balloon 
symmetric about y = 7r/2, while (iii) has flow blocked 
along 7r/2, since the upper half poles are both rotating 
together to 7r/2. This is unphysical and non-generic: the 
rotation mechanism of Sect. 4 can only lead to this under 
extra, nonphysical, symmetry, which of course is exactly 
what the method of Sect.2F creates. With two generic 
poles, case (iii) would not have occurred. Regrettably, 
this generic version is not analytically tractable. 

However, while it turns out that type (iii) never en- 
counters finite time singularities, not all the "good" 
types, (i) and (ii) trajectories are free of disease. That 
is, there is, for each a and 5, a minimum gap, 22/™™ 
between the stagnation points that allows the interface 
to squeeze down through the gap, and then re-emerge, 
blooming out into a balloon. For any smaller yo, the in- 
terface is squeezed into a cusp, unable to pass through 
the gap without penalty of a singularity. This attempt 
to squeeze through and balloon out is the generic disease 
that (7 = theory is plagued by: there are a fraction of 
initial conditions that fail. In fact, this is precisely where 
surface tension needs to be enlisted. With surface ten- 
sion the stagnation points are no longer constants of the 
motion, and indeed will move apart just enough to allow 
the incipient balloon to pass through the gap. It is note- 
worthy that in experimental studies such a phenomenon 
always appears at the initiation of flow (cf. E]j|]). 



So far we have considered a channel filled with fluid 
infinitely far downstream. This is of course unphysi- 
cal. Any experimental apparatus introduces by neces- 
sity some additional boundary condition on the physical 
fluid far downstream, requiring mathematical boundary 
conditions to model this termination. 

We recall that the possibility of adding a sink located 
at some flnite position was discussed in Sect. 2F. We 
can think of other ways to have the fluid itself flnite. 
First, consider the idealized Hele-Shaw cell. At a long 
distance downstream we erect a baffle cross-wise to the 
channel - say at a; = 0. Behind the baffle we have a 
pump controlled to maintain an exactly constant unit 
flux of fluid through the baffle. With a uniform enough 
baffle, we have v{0, y,t) (x x as an approximate boundary 
condition. Thus Vy = = ~dyP on a; = or p{0, y, t) = 
—Pg{t), Pg the positive gauge pressure on the flnite fluid 
from the interface at p = to the baffle. That is 



Re/ = atReC = ^ = C<,(0 • 



But with / reflection symmetric, we have 
.f{C,t) + f{Ct) = atC + C = 2C<,(i) 



(115) 



(116) 



Assuming the fluid is analytic over any region contain- 
ing ReC = in its interior, we then have by analytic 
continuation 



/(C) + f{2^g -0^0, 



(117) 



for all C in the region of analyticity. This exposes the real 
power of reflection symmetry: not only is there a relation 
of the upper physical channel to the lower unphysical 
one, but under flnite boundary conditions from very high 
pressures to very low ones. In this case there is no full 
exponential decoupling of efflux from interface motion. 
This is precisely the "enforced symmetry" between C and 
— C mused about in the sink solution of (|4^) with a = 1 in 
Sect.2F with its upstream pointing 1/2 SaflFman- Taylor 
flnger. We will explore this momentarily, after discussing 
the variant to Hele-Shaw, and a related other pair of 
terminations. 

An obvious variant to flxed velocity on the cross- 
channel line at x = is to simply open (cut off the end of) 
the channel, so that p(0, y, t) = pa =const=atmospheric 
pressure. We then have 



Pa ^ , Pg _ 
V{t) " 

Just as before, we now have 



Re/ = at ReC = = ^ e,(/) > . (118) 



/(C) + /(2G -0 = 



(119) 



so that both variants entail the identical calculations, 
save for the driving fluxes: 
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2 = /'(C)/*(-C) + /'(-C)A(C) 

in the Hele-Shaw case, whereas 



2Vit) 



2Pg 



(120) 



(121) 



in the constant pressure termination, ultimately deter- 
mining the non-steady V{t) in this case, as we saw in 
the most elementary versions of (119) with f = f3 + ( in 
(l)-(ll) of 2F. 

The other pair of variants replace the cross-channel 
line at 2: = with a small circular aperture of radius a 
all along which either Vr = — l/2a so that F = 1 by (0), 
or again p = Pa and Re( = ^g{t). These circular aperture 
problems are mathematically related by exponentiation 
to the cross- channel line versions, and technically much 
harder to discuss with closed solutions. However with 
p = Pa on the circular aperture, there must be a singu- 
larity in the interior of the aperture to sink the full-flux 
that must enter it if we seal off the channel arbitrarily 
far downstream, so that all fluid must efflux through the 
aperture, and in this case, flow stagnates far to the right, 
and so whatever we do far enough to the right will in- 
deed be exponentially suppressed. If the singularity is 
just a simple pole, then it is a sink, generally moving 
within the interior of the aperture. By circle symme- 
try, the analogue for (118) is the fluid gathering into a 
moving sink to right of x = 0, rather than becoming flat 
at infinity. For example, fluid with surface tension after 
emerging from the shaping channel would form a vena 
contracta, and so, reminiscent of a moving sink. This 
was the physical motivation of our consideration of ( ^ ) 
with a = 1 in Sect. 2. F. 



Let now attempt to solve for an / obeying (117) or 



(119). Setting C^C + ^s, M is 



It is easy to check by direct substitution that 

with A arbitrary. Consider 

AiO^^ + Y.^klnil-e'^^-'^) , 

and so, 

/(C) = C - Ut) + E 1^(1 - eC.+««-C) 
-^afcln(l-p>-«''+^) . 



(122) 



(123) 



(124) 



(125) 



Eq.(125) is the entire class (|27| ) of solutions meeting our 
boundary requirements (/3 = —S,g)- 

As a first exa mple, consider just one a and choose Co = 
—(.g + iTT- (125) then is 



/(C) =C-Cg + a ln(l + e-^) - a ln(l + 6'^^^+'^ ) 



For ^g this solution is a single Saffman- Taylor finger 
with an arbitrary width. We insist however that there be 
no flux going off to inflnity, in fact no flux for ReC > 2Cg 
for a fully pinched vena contracta. (This would have been 
automatic in the case of the circular aperture.) To sink 
all flux requires a = 1 (cf. the discussion after Eq.(^5|)), 
and so 



/(C) = C-e5 + ln(l + e-«)-ln(l + e 



). (127) 



which for Cg ^ 1 is precisely a A = 1/2 Saffman- Taylor 
finger (^). This is our first piece of evidence that A = 
1/2 is connected to finiteness. 
But, all is not well. Notice that 



/' 



1 



1 



1-l-e-C 



e "^ 9 



/'(-(») = 0, (128) 



and 



/'(- 



1-1 = 



(129) 



But then, unless ^g is always infinite, (f2Q) and (121) are 



only compatible with V{t) = 0, and so these solutions 
are purely static and not what we seek. 

Consider then more a^. By Eq.(119) if Cfc is a singu- 
larity of /, so too is 2C(, — Cfc- But then, each asymptotic 
stagnation point condition (p3) becomes two conditions: 



/(-Cfc) = ^fc, and /(Cfc - 2C3) = 4 



(130) 



(126) 



Together with the C — > 00 equation for /?, there are about 
twice as many equations as variables unless = 0, in 
which case Pg = 0, and there is no motion. We have 
already seen that one real C has no flux, and it is rea- 
sonable clear that all other cases entailing too many 
equations are inconsistently over-determined. By ( |122| ) 
f'{—oo) = ^ /'(oo) = 0, and so there can never 
be flux with A = 1/2. The above comments of over- 
determination hold for all A. That is, our first four 
schemes of finite termination allow no motion for /'s of 
class ( p7| ) with any finite number of singularities. 

On physical grounds, the fluid emerging into atmo- 
sphere becomes 3-dimensional, and the derivation of 
Darcy's law breaks down. Eq. (118) must be too strin- 
gent. Equivalently, it is not feasible to have a baffle with 
V (X X all along its length. To the contrary, we easily 
imagine fluid racing vastly faster through some holes in 
the baffle rather than others; this choice can readily vary 
in time under minor perturbations of the pump action, 
etc. So, there are hosts of singularities very close to the 
line Rez — 0, and (117) fails for failure of analytic con- 
tinuation. (This is most probably an over-exaggeration: 
it seems not necessary that Rez = is truly a natural 
boundary.) On reflection, these comments imply that the 
physical experiments that have been performed contain 
dynamically determined analyticities, and so are incom- 
pletely posed boundary data configurations. 

Let us consider a fifth scheme of finiteness, of a totally 
different character from the previous four. Consider an 
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infinitely long channel, only partly filled with a finite 
body of fluid, with p = on the left driven face, and 
p = pi = —pg, Pg > on the other, right, free interface. 
The equation of motion for / on the left face, are as usual 
(Eq) with V surely non-zero: 



2V{t) = f{C)M-0 + f'{^C)ft{C) 



The right interface lies at 



ReC = 



Pa 
V{t) 



(131) 



(132) 



and so, with ( C(Coji), Co Lagrangian coordinates, 
with free interface transported to itself, 



ReCt = L 



(133) 



and so by ^) 



V{t) - l/f 4 + Ref'ft on C + C = ■ (134) 

By reflection symmetry, we then obtain a second field 
equation in consequence of the second free interface: 

2v{t) = 24/'(c)/'(2e, - c) + noM2Cg - c) 

+ /'(2^s-C)/t(C) • (135) 



The fluid must now simultaneously obey both pde's, 
(131) and (135). It is unquestionably true that this 
system must have solutions of a physical character, as 
otherwise the entire 2-d theory should have to be dis- 
carded: This fifth version of finiteness is entirely well- 
posed within a conformal 2-d context. Singularity struc- 
ture, of course, is more subtle than our considerations so 
far, but nevertheless if the right interface needs to be a 
natural boundary (i.e. no further analytic continuation 
possible), it is surely the case that so too must be the 
left, because the physics at both are identical. 

Let us now consider a class ( p7| ) solution. (Imagine 
the fluid initially in such a state of perfect repose that 
its / can be naturally analytically continued to+oo.) 



In this case we can take the limit of (131) and (135) as 
ReC oo. We then deduce that 



4/'(+^)/'(-o«) = • 



(136) 



But /'(-f oo) = 1 (dass (|^) and F 7^ ^ 4 ^ 0' and 
so we have class ( p7|) with 



/'(+^)-l, /'(-oo) = 0. 



(137) 



and so only A = 1/2 solutions. 

Suffice it to say that the equations of motion do pos- 
sess a well-behaved solution and we may conclude that in 
the only well-posed 2-D finite system we can construct, 
finiteness alone determines pattern selection. 



7. DISCUSSION 

We have returned to the Saffman- Taylor problem with 
the viewpoint of it as a dynamical system in order to 
better understand the evolution of its solutions. In do- 
ing so we have carefully re-thought the relevant bound- 
ary geometry and conditions and realized that reflection 
symmetry rather than periodicity is to be imposed. This 
led to two significant consequences. 

First, reflection symmetry and analytic continuation 
naturally promoted the equations of motion from a re- 
lation pertaining purely to the interface, to one of a 
field character throughout the fluid. In consequence, we 
need never consider the usual Hilbert transform bound- 
ary methods, instead directly, and largely algebraically, 
obtaining solutions and their dynamics. 

Secondly, the fluid equations naturally link /(C) and 
/(— C), so that very far downstream details of termina- 
tion potentially couple to the very far upstream (above 
physical fluid) details, such as the singularities determin- 
ing the flow. This counter-intuitive failure of termination 
details to exponenitally decouple from the behavior of the 
interface propelled us to contemplate that finiteness in 
this problem is apt to be a deeply significant "singular" 
perturbation upon the "physics" of the infinite channel 
problem. In consequence, we formulated the theory from 
the beginning to include the possibility of variable flux, 
a necessity of finite configurations. 

Employing our reflection-symmetric field equations, we 
readily produced a variety of elementary solutions and 
then the pole-dynamics family (|27|). In particular we 
determined the general form of all translation-invariant 
solutions, and the simplest pole- type solutions more com- 
plicated than the original Saffman- Taylor class. These 
are characterized by a downstream sink, which when fully 
sinking all flux, produces an upstream pointing 1/2 finger 
surrounding the zone of efflux. Considering how the na- 
ture of this flnger is contingent upon final termination (a 
full sink at finite distance), and considering the ±C sym- 
metry of the equations of motion, it is impossible not to 
wonder that we might be touching upon the origin of pat- 
tern selection. In this context the reader should not be 
troubled by the circumstance that the sink is now within 
the body of physical fluid. He should not (or should) be- 
cause this is identical to the situation in the usual infinite 
channel flow, when the full sink instead of appearing at 
finite ReCs is at ReCs = +00, still fully within physical 
fiuid. (One might think of a Mobius transformation ro- 
tating the point at infinity to a proximate point.) This, 
in fact, should trouble the reader, because it means that 
efflux in the usual case (ReCs = -l-oo) has not been phys- 
ically treated: The volume of fluid is conserved only be- 
cause 00 — 1 = 00. As we shall see in [||, a full treatment 
of all real fluid has significant physical consequences. In 
particular, it transpires that each unstable mode /„ of 
( ^ ) requires power from the energetic sources driving 
the flow, so that under pump control, the exponentially 
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growing modes are sharply supressed, leaving behind, at 
best, resummations such as class (^7|). 

We proceeded to analyze the evolution of an arbitray 
flow, although largely within the context of class (^7|), to 
better understand how wcU-formualted the theory is, and 
some general boundary violating circumstances of finite- 
time singularities - namely those that have been put in 
evidence in the prior literature. We later went on to ex- 
hibit the general circumstance of a finite-time singularity 
within class (p7|), which is the situation of an incipient 
balloon attempting to negotiate passage through a pinch- 
ing pair of stangation points. With arbitrarily small sur- 
face tension, the class ( p7| ) flow is unaltered until the 
tip of the penetrating fluid is approaching a cusp with 
diverging curvature. At this point the singular pertuba- 
tion renders the stagnation points no longer constants of 
the motion. However, as soon as the pair has separated 
far enough to allow the balloon to form, the curvature is 
quite finite, and class (|2^) is again correct, save that the 
stagnation points are just far enough apart to allow the 
minimum waisted balloon to pass: Had we chosen ini- 
tial data to have been these new locations of stagnation 
points, the (7 = theory would have fully sufficed, and 
produced the physical solution. Following the last para- 
graph of Sect. 5, taking the simplest class ( |27| ) solution 
with one pair of complex C^'s with Recife = 1/2, one can 
flnd that Ima/c with the narrowest waist, and observe the 
strong similarity between the asymptotic balloon and the 
best developed experimental one of jl]. 

We next observed, purely within class ( ^ ) however, 
that with singularities coming in clusters, each cluster 
well-separated in ReC from another, that the solution has 
asymptotic regimes, with singularities far to the left play- 
ing an exponentially insignificant role upon the shape 
of the interface, while all the others are very close to 
ReC = 0, and as we demonstated, having migrated to 
lm( = or TT. That is, until another cluster of singular- 
ities at the left arrives close to ReC = 0, at which time 
it joints into the asymptotics of those already there, the 
interface evolves as a single Saffman- Taylor finger. As 
another cluster arrives, that finger metamorphoses into 
another of a new width if y^' a k of those arriving dif- 
fers from zero. That is, class (^7|) has the asymptotics of 
always a single finger, but of generally metamorphosing 
width. To establish A = 1/2 on these grounds requires 
a reason for afe = 1 for just those singularitieis near 
ReC = during the period of time of the physical ex- 
periment: with ^ afc = 1 for all singularities, including 
those arbitrarily far to the left, demonstrates nothing 
about physical pattern selection. What has fundamen- 
tally characterized Sect. 3-5 is our focus upon temporal 
evolutions, accomplished via class (p7|), by considering 
this flow as a dynamical system. 

Finally, we pick up on the ±C symmetry, a consequence 
of dynamics imposed upon ReC = of a reflection sym- 
metric system. With any boundary fixing on another 
curve, say ReC = ig{t), a sharp relation of -|-C to — C must 
follow, such as with fixed pressure along a downstream 



line, yielding (117). This makes it clear that arbitrarily 
far downstream terminations {^g ^ 1) somehow become 
entangled with ReC ^ —1, the domain of singula ritie s 
that determine the shape of the interface. Although (117) 
allows of no class (p7|)^ solutions, meaning including 
singularities far to the right as well as those to the left 
of ReC = 0, this does not mean that theere are no solu- 
tons: this is largely the insufficiency of finite order class 
(p7|), even when extended to include higher order singu- 
lar terms. (We shall see this in 1^.) However, it is clear 
on physical grounds that (117) is too stringent a symme- 
try, and is to be replaced by myriad singularities in the 
flows's analytic continuation beyond termination. This is 
a serious modiflcation of the problem, since these singu- 
larities are dynamical and of a priori unknown character 
and locations, instead determined by all the mechanical 
vagaries of the innards of a pump and so forth, and so no 
longer a physically sensibly posed problem. Instead, the 
physical problem is one of boundary geometry over just 
the experimentally observed body of fluid, and hence is 
one of incompletely posed geometry and data. This en- 
tails solutions no longer unique, requesting a physical 
mechanism to select among branches etc. 

Accordingly, within the machinery and formulation at 
hand, just one choice lay open, which is to consider 
the purely 2-D conformally well-posed problem with two 
free interfaces. The full treatment of this purely non- 
autonomous, non-periodic problem is the subject of 
We reflected some of its introductory matter into the last 
paragraphs of this paper to complete the flow of our con- 
siderations. It is worth mentioning that there is a class 
(^ ) solution with precisely on real Cfc (and ak = 1), and 
no others whatsoever within class (p7|). 
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